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Abstract 

We investigate the supercooling of a nematic liquid crystal using fluctuat- 
ing non-linear hydro dynamic equations. The Martin-Siggia-Rose formalism 
is used to calculate renormalized transport coefficients to one-loop order. 
Similar theories for isotropic liquids have shown substantial increases of the 
viscosities as the liquid is supercooled or compressed due to feedback from 
the density fluctuations which are freezing. We find similar results here for 
the longitudinal and various shear viscosities of the nematic. However, the 
two viscosities associated with the nematic director motion do not grow in 
any dramatic way; i.e. there is no apparent freezing of the director modes 
within this hydrodynamic formalism. Instead a glassy state of the nematic 
may arise from a "random anisotropy" coupling of the director to the frozen 
density. 

1 Introduction 

The study of a supercooled nematic liquid crystal and the possible formation 
of a nematic glass^ is potentially richer than corresponding studies of super- 
cooled simple fiuids. The presence of anisotropy in the nematic liquid due 
to the overall alignment of the molecular long axes introduces orientational 
degrees of freedom into the description of the system and yields a model sim- 
ilar to a spin glass with translational degrees of freedom coupled to the spin 
fiuctuations. We use the term "nematic glass" to describe a liquid crystal 
where both the translational (density) and orientational (director) fiuctua- 
tions are frozen. Assuming that we supercool the liquid starting from its 
nematic phase (rather than the isotropic phase) we expect a nematic glass to 
have long-range orientational order and thus be similar to a mixed magnetic 
phase where both ferromagnetic and spin glass order coexist.^ The presence 
of both translational (density) and orientational (director) fiuctuations and 
their coupling leads to a novel glass-forming system.^ Glass formation in this 
system could potentially occur in a two-stage process where the density or 
director modes freeze first, followed by the other, or in a process where both 
freeze simultaneously. 

In recent years a theoretical approach to the study of glass formation 
has been developed using either mode-coupling calculations'^ or fiuctuating 



nonlinear hydrodynamics.^ These approaches were initiated by Leutheusser^ 
who showed that a model of a dense fluid obtained from kinetic theory ex- 
hibits a sharp glass transition where the system becomes nonergodic. Das, 
Mazenko, Ramaswamy and Toner^ developed an equivalent model on the 
basis of fluctuating nonlinear dynamics. Subsequently, Das and Mazenko^ 
discovered a nonhydrodynamic mechanism which cuts off the sharp transi- 
tion, leading to a rounded transition. These authors claim that the cutoff is 
due to the proper mathematical treatment of the relationship P = pV where 
P is the momentum density, p is the mass density, and V is the velocity field. 
The physical origin of this cutoff is unknown. More recently, Schmitz et al.^ 
have argued that the calculation of Das and Mazenko is not correct. Fur- 
thermore, they claim that the perturbative calculations of refs. 5 and 6 do 
not properly account for detailed balance. Restoring detailed balance to the 
perturbation theory apparently restores ergodicity and leads to a rounded 
transition. However, Das and Mazenko and Schmitz et al. all agree that 
the original Leutheusser theory as well as its hydrodynamic version with no 
cutoff mechanism are a good approximation to the more complete theories 
of refs. 6 and 7 in describing the growth of the viscosities in the pre-glass 
transition regime. Eventually the growth of the viscosities is limited by one 
of these cutoff mechanisms. Comparison of these theories with experiments 
yields some encouraging agreement,^ though this agreement is by no means 
complete, especially at very low frequencies. 

In this paper we study the formation of a nematic glass using fluctuat- 
ing nonlinear hydrodynamics. We will limit our attention primarily to the 
pre-glass transition regime in light of the discussion above, and we will not 
concern ourselves with the question of detailed balance or the Das Mazenko 
cutoff mechanism. The primary advantage of the fluctuating hydrodynamic 
theory is that new slow variables (such as the director modes or variables 
associated with broken translational symmetries^) are readily incorporated. 
As described in detail in the next section, we supplement the nonlinear hy- 
drodynamics equations used by Das et al. to study simple fluids, with two 
equations describing the dynamics of the director in a compressible nematic. 
The full set of equations we employ also includes a nonlinear coupling {n-\jpY 
between the density and director n. Density fluctuations couple to some, but 
not all of the nematic viscosities through this coupling. The Leutheusser 
feedback mechanism then leads to the enhancement of these viscosities with 
a universal prediction for their power law behavior. However, there is appar- 



ently no director feedback mechanism within this formahsm, and we find no 
evidence for the freezing of the director modes. Nevertheless, it is tempting 
to speculate that the term {n ■ \j5pY in the free energy leads to freezing of 
the director. Once the density has frozen, or almost frozen this term will 
mimic a random anisotropy field in an amorphous magnet. ^° The frozen den- 
sity gradients \jp play the role of the quenched random axis. This model is 
believed to exhibit a spin- glass phase, which would correspond in our case to 
frozen director modes, with no nematic long range order in an infinite system. 
However, for a finite sized system, there would be apparent long-range order, 
especially if the coefficient of the biquadratic term is small. Whether this 
proposed freezing of the director takes place immediately upon freezing of 
the density or requires further supercooling is unclear and beyond the scope 
of our present theoretical treatment. 

This paper is organized as follows: in the next section we formulate the 
nonlinear hydrodynamic equations for a compressible nematic including the 
above-mentioned nonlinear coupling of density and director modes. In sec- 
tion III we use the Martin-Siggia-Rose formalism to study the effects of the 
nonlinearities on the bare transport coefficients. Finally in section IV we dis- 
cuss the implications of our calculations for the growth of the viscosities and 
the mode structure as the nematic is supercooled. Various technical details 
appear in the appendix. 

2 Non-Linear Hydrodynamic Equations 

In contrast to the spherical molecules of simple liquids, the molecules of liq- 
uid crystals are elongated in shape. ^^ Intermolecular interactions cause the 
anisotropic molecules to align along a preferred direction denoted by the 
director vector n(x, t). Fluctuations in the director can extend over macro- 
scopic distances and decay over ffnite times, leading to new hydrodynamic 
modes in addition to the shear and sound wave modes of a simple liquid. 
The broken rotational symmetry also permits more viscosities: in a com- 
pressible nematic liquid crystal, there are six independent viscosities.^^ The 
linear equations for the hydrodynamic modes of a liquid crystal have been 
known for almost twenty years. ^^ A systematic method for deriving the non- 
linear contributions is to write them in the form of generalized Langevin 



equations: ""^"^ 

di,./dt = v.m-Y.J d'x' r,,(x') ■ j^^^ + 0„ (2.1) 

where %l)i{^^t) represents one of the seven possible hydrodynamic fields: the 
mass density p(x, t), three components of the momentum density P(x, t), 
the energy density e(x, t) and two components of the fixed-length director 
n(x, t). The label i in this equation denotes the type of field, as well as the 
vector index on P and n. Vi [ip] represents the reversible part of the dynamic 
equations and is given by: 

n^i^)] =T.J d'x'iM^), U^')}j^i;^y (2.2) 

and H is the energy obtained by integrating the free energy density F[iIj]: 

H = I d^xF[%l){^)]. (2.3) 

The Poisson bracket in eqn. (2.2) is defined in its usual manner as: 

5ipi 6ipj 5ipi 5ip 

a,l3,k ^^fc ^-^fc '^^k '^^k 






where r^ is the k^^ component of the vector r"^ which points to the (3*^ atom 
of the a*'* molecule of the liquid crystal. The second term on the right hand 
side of the equation of motion (2.1) represents the dissipative contributions in 
form of the dissipative matrix ry(x, x'). Finally, Oi(x, t) denotes a Gaussian 
noise source which satisfies: 

(Gi(x,t)e,(x',t')) = 2A;BTr,,(x-x')5(t-t')- (2.5) 

It has been argued that the dominant transport anomalies at the glass tran- 
sition in a simple liquid are due to the slow decay of density fiuctuations, 
and the effects of energy fluctuations can be ignored. ^'^ We shall denote the 
global preferred direction of orientation by no = e^. Local fluctuations in 
orientation are specifled by 5n = n — no where we restrict 5n = {drix-, Suy, 0) 
to linear order in Sui}^ Then n satisfles n- n = 1 up to order ((5n)^. Further- 
more one can assume that the wave vector k of a disturbance lies in the x — z 



plane. This allows us to treat drix and 5ny as longitudinal and transverse 
fluctuations with respect to k respectively.^^ Our set of dynamical variables 
%lji then includes the density, three components of the momentum, and the 
director fluctuations, n^ and Uy, which are respectively equal to 5n^, and 5ny 
in the present approximation. 

In order to evaluate the Poisson bracket between the six dynamical vari- 
ables p, Px,Py,Pz, Ux and Uy we need a microscopic description for each vari- 
able. We define: 

p(x,t)=^m^5(x-r"^(t)), (2.6) 

P,(x,t) = ^i^"'^(t)5(x-r"/^(t)), (2.7) 

n,(x,t) = ^^<(t)5.,R.. (2.8) 

We model an elongated molecule with two atoms only: R" is the center 
of mass of the a*'^ nematic molecule and n" = r^ / | r° |, and r" is 
the relative vector between the atoms, which, we shall assume, model an 
elongated molecule of mass m°. The variable (3 in r"'^ takes on values 1 or 
2, and r" = r"^ - r"^ while R" = (r^^ + r"2)/2. 

With the definitions (2.5) - (2.8), the Poisson brackets required in (2.2) 
and defined in (2.4) can be computed in a straightforward manner with the 
following results: 

{^(x), P,(x')} = VM^)K^ - X )); (2.9) 

{P,(x), P,(x')} = -V,(5(x - x')P.(x)) + V:(5(x - x')P,(x')); (2.10) 

{P,(x),n,(x')} = [(A + l)5,,nfc(x')/2 + (A - l)4jni(x')/2 - A(n,n,nfc)(x')] ■ 
Vfc[5(x - x')] + 5(x - x')(V:n,(x')). (2.11) 

where Vj = d/dxi and V^ = d/dx'^ etc. In writing eqn. (2.10) we have 
introduced the "form factor" A which is related to the shape of the molecules 
and equals unity only in the limit of infinitesimally thin molecules. ^^ All other 
Poisson brackets are zero. 



The final step in evaluating Vilip] is to calculate SH/Sipj. Expressing the 
free energy density F{p, P, n) as a sum of kinetic and potential energies, we 
rewrite (2.3) as: 

H = jd'x{e, + ei + el + el) (2.12) 

where Sk is the kinetic energy density of the molecules, e(^ is the potential 
energy density due to density fluctuations and e" is that due to the director 
fluctuations. Finally e^ is an energy density due to the coupling of density 
and director fluctuations. For the kinetic energy density we have: 

efcx,t=y^, 2.13 

2p(x,t) 

For the potential energy density of the density fluctuations we choose the 
simplest form incorporating \/p: 

5S(f , t) = ^{6p{x, t)f + |(Vp(f , t)f (2.14) 

where A and B are phenomenological constants, and 5p = p{x, t) — po with 
Po being the uniform density. The gradient term in (2.14) is rotationally 
isotropic which is unrealistic in an anisotropic system like the nematic. We 
will incorporate the effects of anisotropy below in e^, the coupling of Vp to 
h 

The simplest choice for e^ is: 

e:i^,t) = iir(V.n,(x,t))(V.n,(x,t)) (2.15) 

where i^ is a Frank elastic constant and repeated indices are summed over. 
In general, the symmetry of a nematic allows for three independent elastic 
constants^^ corresponding to the distortions splay, twist and bend. We have 
made the simplification of setting the three elastic constants equal to K. 
While this equality is broken upon renormalization (see Section III), the 
difference between the elastic constants is not large and we will ignore it. 

Finally in e^ we incorporate the expected anisotropy in density fluctua- 
tions and choose, 

<(x,t) = i/(n(x,t)■Vp(x,t))^ (2.16) 



where / is a phenomenological coupling constant. In general one might in- 
troduce a couphng of the form |(n ■ Vp)^ + -^(n x Vp)^ to account for the 
expected anisotropy in the density fluctuations, i.e. the energy of such fluc- 
tuations should depend on the relative orientation of the director and the 
wavevector of the fluctuation. However, because of the following vector iden- 
tity: 

(Vp)2 = {n ■ Vpf + (n X Vpf (2.17) 

we can eliminate the coupling proportional to /' in favor of a redefinition 
of / and the inclusion of the last term in (2.14). This choice simplifies our 
perturbation theory in section III, and allows / to be negative. However 
rodlike molecules will probably be characterized by positive values of / since 
the director will prefer to align perpendicular to the wavevector of the density 
field. We can think of (2.14) and (2.16) together as providing a simple model 
for the static structure factor of the nematic, where density correlations are 
not isotropic in space but depend on the local director orientation. The 
effect of local structure on a simple fluid as it is supercooled was considered 
by Das.^^ 

Using (2.9) - (2.17) we can evaluate the reversible terms in (2.1). For the 
P equation of motion we can express the result conveniently in terms of the 
divergence of a reactive stress tensor, i.e. 

V>. = - ^ VjctJ^, i,j = x,y,z (2.18) 

j 

where 



P,R 



+ p/Vfc(nfc(n-Vp)) + |(n-Vp)2" 

+ I{n-Vp)njVip + K{Vine){Vjne) 

- Efe Prjk{-KV^nk + J{n ■ Vp)Vkp) (2.19) 



and 



1 1 

Pijk = -(A + l)njSik + -(A - l)niSjk - XriiUjUk (2.20) 

The stress tensor appearing in (2.18) is not symmetric. However, the equation 
of motion for P is not sensitive to this asymmetry. A symmetric choice for a^ 
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would yield the same forces as well as guaranteeing conservation of angular 
momentum^^ which is not an issue in our analysis. 

The reversible part of the director equation of motion is given by, 

Vn^ = — P ■ Vrii + [Vty. n)i + \ ^ {ukSij - njnjnk)Akj i = x,y (2.21) 

" k=x,y,z 

where 

Akj = ^{VkPj+VjPk) (2.22) 

Zp 
and, 

n = —{y X P) (2.23) 

A dissipative stress tensor can be introduced to write the dissipative con- 
bution to (2.1 
via the relation. 



tribution to (2.1). The dissipative momentum stress tensor a^j is defined 



j=x,y^z^ j=x,y^z J 

P D 

The uniaxial symmetry of the nematic dictates that ctj^' has the following 
form,^^ 



a. 



ij 



2(z/° + z/2 - 2u^)nknmAmkninj + 2{u^ - U2){ninkAkj + riinkAik) 

+ {l^5 -^1 + V2){'^i^j^kk + ^ijUkUmAkm) + ^V^Aij + (z/° - Z/°) Afefc5i/2.25) 



The five bare viscosity coefficients z/°, z/^ , z/g, z/^, and z/^ are in general not 
equal to each other. If z/° = z/2 = z/g = z/4 — z/5 we recover the usual form for 
the dissipative stress tensor of a simple fluid. If the nematic is incompressible 
then all terms in a^j which are proportional to Akk mush vanish, implying 
that z/2 = z/4 and z/5 = 0, and three independent viscosities remain. As we are 
interested in density fluctuations we will work with the full tensor displayed 
in (2.25). Comparing (2.24) and (2.25) we can identify the nonzero elements 
of the viscosity matrix Tpipf 

^p.p. = -{< + 4)^1 - ^sV^, (2.26) 

^PyPy = -4^1 - ^3 VL (2.27) 



TP.P. = -i^lVl - (2Z.? + z.° - z.° + 2z.°)V^, (2.28) 

^p.p. = -K + ^3°)V.V.. (2.29) 

All other elements of the viscosity matrix TpiPj = 0. If we go the limit 

of a simple fluid (z/J* = 1^2 — ^3 — ^4 ~ ^5) ' then TpiPj reduces to the tensor 

Lij, used in refs. 5 and 18. 

The dissipative contribution to the director equation is proportional to the 
"molecular field", 6H/6ni, with the proportionality constant conventionally 
written as I/7?, where 7^* has the units of viscosity. Thus we identify, 

Tmnj = —Sij, ij = x,y (2.30) 

7i 

All other elements of the viscosity matrix T^ are zero. In particular all 
off-diagonal elements r„-p. can be shown to be zero on the basis of time- 
reversal symmetry. ^^ 

The compressible nematic is thus characterized by six hydrodynamic fields 
and six independent viscosities: i''l,i>2,i'^,i'4,i'^ and 7°. Another viscosity 
72 commonly discussed in the literature which describes the torque exerted 
on the director by a shear flow is related to 71 via 72 = —X'^i where A is the 
reactive coefficient introduced in eqn. (2.11). 

Our final nonlinear hydrodynamic equations for the compressible nematic 
are as follows: 

^ + V-P = (2.31) 

dP 

^ = - E V,4' ^ = x,y,z (2.32) 

3ti ■ -> 1 -. -. 

ot p 

+X Yl Akenk{6ii - niUi) 

k,£=x,y,z 

K I ^ 

-\ — frV^ni fr(n- Vp)Vi/), i = x,y (2.33) 

7i 7i 



where 

with a^j and a,ij given by (2.19) and (2.25) respectively. 



P P,R I P,D /o Q/l^ 



9 



3 Transport coefficients 

To investigate the effects of the nonhnearities in the equations of motion 
(2.31) - (2.34) on the transport properties of the nematic we use the Martin- 
Siggia-Rose (MSR) formahsm.^^ This formahsm allows us to calculate the 
correlation functions Gij{x,t;x',t') =< 5ipi{x^t)5ipj{x' ,t) > where ipi repre- 
sents any of our six hydrodynamic fields and the brackets refer to the average 
over the noise source Oj defined in (2.5). It will also enable us the calculate 
response functions, and obtain corrections to the bare viscosities introduced 
in the previous section. We refer the reader to refs. 6 and 19 for full details 
on the MSR method; we summarize the essential ideas here. 
We define a generating function Zu[i/j,ip] as follows: 



(3.1) 
where C is a constant and ip collectively represents the six hydrodynamic 
fields ipi. Functional differentiation of Zu with respect to U generates the 



correlation functions Gij: 



6HnZ^ 



u 



^''^'^'^''^''^'^~ 6U,ix,t)6U,ix',t') ^^-^^ 

The six auxiliary fields ip were introduced to exponentiate each of the six 
hydrodynamic equations of motion (2.1). The integration over the noise 
source Gj has been replaced by an integration over the fields ipi. The action 
Au is given by: 

A(/ [■?/',■?/;] = dxdtdx'dttpi{:>c,t)(3~^Tij{-x.,t;x',t')t[jj{'x';t) 

+t f d^d4,{^,t)[^^^jl^ - H,[^]]. (3.3) 

where Hj is defined by: 

^.M = V,m - E/^xWr,,(xY)^^^-^ (3.4) 

Introducing the "vector" = {p,P,nx,ny,p,P,hx,hy), we can then rewrite 
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Ajy [■?/',■?/'] = Au[(p] as follows: 



i I dld2 ^ M^)Go\l, 2)0^(2) (3.5) 

+ E If [ dld2...dN Y. Kbi 0.(1)0/3(2)0,(3).... 

a/57 . . . ^^^ 1^ ' 



Here the integration variables 1, 2, 3, . . . stand for (x, t). The explicit expres- 
sions Gq^ for the vertices Va/?^... will be presented later. The action Au gener- 
ates the full non-linear hydrodynamic equations for a nematic liquid. When 
we omit the vertices what remains then is a pure quadratic Gaussian field 
theory. The corresponding Ajj contains the inverse of the linearized correla- 
tion matrix Gq^ and thus generates the linearized hydrodynamic equations 
of a nematic liquid. We first discuss this limit before doing perturbation 
theory in the vertices. 

The linearized theory, as well as the perturbation theory, are most easily 
discussed in Fourier space. We Fourier transform (3.5) in space and time 
recalling our discussion following (2.5) where we assumed that the wavevector 
k lies in the x — y plane. The matrix GQ^{kuj) is displayed in Table I where 
we have introduced the following tensors for notational convenience: 



K 

Y 

K 



^ij = ^ij — {^+'^)kz 



+SixSjzYi^ ~ ^)^^ (3-6) 



+ (5,,(5,,rL + {5,x5,. + 5.,5,,)r°, (3.7) 

where F^^, V^y, F°^ and F°^ are the Fourier transforms of the viscosity matrix 
elements (2.26) - (2.29) respectively. 

The inversion of Gg is most easily accomplished by decomposing P and 
n into longitudinal and transverse components. The calculation is straight- 
forward though tedious especially for the longitudinal portion. The elements 
of the matrix Gq provide the physical correlation and response functions for 
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the linearized theory. Correlation functions of any two hatted variables are 
identically zero. Correlation functions of unhatted variables can be found 
to leading order in /c^ from the corresponding response functions via the 
fluctuation-dissipation theorems: 

< pp >= -2/3"^x-^m < pp> (3.8) 

< n:j:n^ >= -2(3~^XnIm < n^h^ > (3.9) 

< UyUy >= — 2/3~^Xn-^"^ < ''^y'f^y > (3.10) 

where x is the static density structure factor given by {A + Bq^)~^ and Xn 
is the static director structure factor given by {Kk'^)~^. We could also write 
corresponding relations for the momentum correlation functions, however, 
they are not needed in elucidating the mode structure and calculating the 
corrections to the transport coefficients. 

The density and director response functions in the linearized theory are 
given by the following expressions: 

< s,Ct.^m^K -.) >»^ ^^^:^f^l^,^ (3^11) 

< Sn,it.mA^t, -.) >"^ ^^^^^;;^l,^,^ (3.12) 



< Sn,(k, ^)in,(-k, -^) >"= ,,. , ,„;:„■ ^.,,,^„. (3-13) 



The results are correct up to terms of relative order k'^. The viscosities 
appearing in (3.11)-(3.13) are given by: 



To 








+ (2z/° + 2z/° + z/° - z/°) cos^ e 


(3.14) 


^1 = 


-- z/°cos^(2^) + (z/° + u°) sm'^{29)/2 


(3.15) 




4 = z/° sin^ e + u^ cos^ e 


(3.16) 
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2\ 




.K 



AKul 



'l + Acos2^)2— + 

Po Poll 



(3.17) 



sj 




2\ 




(A 



Po Poll , 



(3.18) 



where 6 is the angle between k and the z axis. 

In writing (3.17) and (3.18) we have assumed that the director modes 
appearing in (3.12) and (3.13) are diffusive rather than propagating. This 
is true for equihbrated nematics^^ where the orientational relaxation time 
of the director (Kk^/'ji)' 
{vlk^/po)-^ and (z/^A;Vpo 
remain true when the nematic is supercooled or compressed rapidly, and 
propagating shear modes can appear. In the absence of supercooling (3.17) 
and (3.18) can be approximated as follows: 



^ is small compared to the shear diffusion times 
r^. We shall see subsequently that this does not 



^ f 



r: 



■po 



^llPo 

VtIPo 



(3.19) 

(3.20) 
(3.21) 



Thus there are two "slow" modes with viscosities r° and r° corresponding 
to the slow relaxation of director fluctuations, while the fast modes with 
viscosities T f and T f are like ordinary shear waves. In addition to these four 
modes we also have two sound modes appearing in (3.11) with speed Cq and 
damping Fq. 

We now calculate the corrections to the linearized theory due to the ver- 
tius V^ appearing in (3.5). If we perform the rescalings ■?/' — * /3"'"^/^-?/' and 
iIj -^ (3~^/'^il) we see that the quadratic part of Au is 0((/cbT)°), and the 
higher-order terms proportional to V^ are of order (kBT)^^^"^^"^ . Thus we 
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can systematically compute corrections to the linearized theory in powers of 
ksT. In particular one-loop diagrams will be 0{kBT). 

The inverse of the correlation matrix for the complete nonlinear theory 
satisfies the formal equation: 

G-\l,2)=G-\l,2)-^{l,2) (3.22) 

which defines the self-energy S. We can then write corresponding equations 
for the renormalized transport coefficients by referring to Table I. As shown 
in ref. 6 the renormalized viscosities are most readily obtained by looking 
at the renormalization of the terms in the action which are quadratic in the 
hatted fields. Thus the elements of the viscosity matrix r^ renormalize as 
follows: 

Tt = K. + ^^P^pSl^) (3.23) 



Tyyik,u;) = T;^ + ^^P^P^{k,u;) 




(3.24) 


TUk,uj) = T:, + f^i:p^p^{k,u;) 




(3.25) 


T,4k,u;) = T:, + ^i:p^pfi,io) 




(3.26) 


lhiik,uj) = ^^ + 


2 7i 2 


~^y{k,Uj) 


(3.27) 


Referring to Table I we find that K/'ji, A and C renormalize as 


follows: 


K 

7i 


(K\ i 

\lioJ k^ 






= 


(K\ i 




(3.28) 


A = 


2 

1^2 




(3.29) 



c' = cl = +^Sp^^ (3.30) 
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The uniaxial symmetry of the nematic as well as conservation of momentum 
allows us to the write momenta self-energies as follows, (cF. (2.26)-(2.29)): 



Ep^p^{k,uj) = - kl{-f2{k,uj)+-f4{k,uj)) 
- kl-f3{k,uj) 

T^ppik^u) = -kl-f2{k,u) - kl-f3{k,u) 



(3.31) 
(3.32) 



Sp^P^(^,^) = -kl-fs{k,uj) 

-kl [27i(A?,c^) + 72(A?,cu) 
-74(A?,c^) + 275(A;,tu) 

Ep p^{k,ij) = -k^k^{j3{k,u) +75(A;,cj)) 



(3.33) 
(3.34) 



where the functions 'ji{k,u),i = 1,...,5 renormalize the viscosities z/j as 
follows: 

Ui{k, u) = u° + -'ji{k, a;), i = 1, . . . , 5 (3.35) 

In the hydrodynamic limit the response functions will have the same form 
as in the linearized theory with the bare transport coefficients replaced by 
their renormalized values which depend on k and uj. As discussed in the In- 
troduction we are ignoring any possible nonhydro dynamic terms which might 
cutoff a sharp transition and focus instead on the behavior of the renormal- 
ized viscosities. Eqns. (3.14) - (3.18) will then be valid for the renormalized 
quantities and we have the following relations for the generalized transport 
coefficients appearing in the response functions: 



r(A;,cc;) = r- p 



^^p^pSkM + ^^P^pfi^^) 



+2^S 



p.p. 



{k,uj) 



(3.36) 



VL{k,^) = ^l- /^ 



^^P.P.(^'^) 
^2 

+;fcf^p.A(^'^ 






^p^pAk,uj) 



(3.37) 
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UT{k,u;) = 4 - /? (1) ^p^p^{k,u;) (3.38) 

Using (3.32) - (3.35) we see that these renormahzed viscosities are well be- 
haved in the limit A; ^ 0. 

The preceding discussion indicates that we need to calculate the follow- 
ing self-energies in order to obtain the renormahzed transport coefficients: 
Sp^p^,Ep^p^,Ep^p^,Ep^p^,Sft^ft^,S„^ft^,Sft^P^, and Sp^^. We have done so 
to one-loop order under the following condition: we keep only those dia- 
grams where there is a possibility of feedback from either density or director 
fluctuations. Thus we consider diagrams where the propagators are either 
Gpp{k,u) or Gninj{k,u),i,j = x,y. (In the Appendix we show that mixed 
propagators cannot yield a feedback). The diagrams will be bubble- type and 
contain two three-point vertices. The symmetrized three-point vertices in 
the action (3.5) that contribute to the diagrams within our approximation 
are given by: 

Va/s-r = - ^1/37(1, 2, 3) + V0Q^(2, 1, 3) 

+V;^,(3, 2,1) + 14^^(1,3,2) 

+V>^,(2,3,1) + V;,^(3,1,2,)1 (3.39) 



where 



\4^,(1,2,3) = ^V;«^(1,2,3) (3.40) 



and the V^'^^ are given in Fourier space by. 



[A6ij{k2 + k^)j + B{{k2 + k3)jk3ik3j 

-{k2 + k^),k2jk^j)]5{l + 2 + 3)(27r)^ (3.41) 

V^^l^{l, 2, 3) = iKY,5^p.5^nAnXi)^k2jkl5{l + 2 + 3)(27r)^ (3.42) 

3I 



V;(5_^(l,2,3) = -(A + l)KY,5^p.5pn,5,n,{i)\k2i + k^i)kl5{l + 2 + 3)(27r) 



(3.43) 
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yijS,^(l,2,3) = -(A - l)KY.5^^^5pn,5,nXi)\k2i + hi)kl5{l + 2 + 3)(27r 



(3.44) 
Fi^j,^(l,2,3) = -i\KY.5^pJpn,5^n,{i)\k2,+h;)kl5{l+2+2.){2T,Y (3.45) 

Fi5,,(l,2,3) = -^J^5„^.5;3Ap(0'(fc3,fc3'J5(l + 2 + 3)(27r)^ (3.46) 
i 

\/i3,^(l,2,3) = -2^(A + l)^5«p/^Ap«'(fc2,A;3.)(A;2.+A;3.)'5(l+2+3)(27r)^ 



(3.47) 



VJJ, ri,2,3) = -z-(A-l)5:<5„,,(5^p(5,p(^)3(%A;3.)(% + A;3,)(5(l+2+3)(2 



TT 



(3.48) 

V^^l^il, 2, 3) = -i/p° ^ (5,p. V^(5^p(i)3(A;2j + A;3,) ■ ((/c2. + A;3x)/c3. + 

i 

(A;2. + A;3.)A;3.)5(1 + 2 + 3)(27r)^(3.49) 

V;^^^!' 2. 3) = i\Uap.5pp5^p{{)\k2, + A;3,)A;2,A;3,5(1 + 2 + 3)(27r)^ (3.50) 



7i p 



Jjp[i) 



k2ik3J{l + 2 + 3){2n)^ (3.51) 

where 6{1 + 2 + 3) = 6{ki + A;2 + k^). The three-point vertices appearing 
in (3.42) - (3.51) arise from the reactive momentum stress tensor (2.19) 
excluding the convective term (PiPj/p). The vertex (3.52) arises from the 
dissipative part of the director equation; as we shall see it ultimately plays 
no role in the growth of the viscosities. 

Using these vertices we have evaluated the relevant self-energies. The re- 
sults are quite complex and are tabulated in the Appendix. The renormalized 
transport coefficients have the following form: 

R r°o r rlP 



T{k,u) = T, + ^P dte'^' j 



2 Jo J (27r)3 

X \Gppik',t)Gppik-k',t)ix-' + fo{k, k') 
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+Gn^n^ {k', t)Gn^n^ {k - k' , t)K^ h{kk' 
+Gn,ny{k\ t)Gn,ny{k - k't)K^ h(k, k' 

+G,,ik',t)Gn^nAk-k',t)PUk,k') 



(3.52) 



UL{k,Uj) 



v°r + - I dte 



Sijjt 



dk' 



2Jo ' ' J (27r)3 
G,,{k',t)G,,{k-k',t)K^g,{k,k') 

Gn.. 0', t)Gn,n, {k-k', t)K^g2{k, k') 



(3.53) 



PT{k,u) 



z/° + ^K^e 



2' 

dk' 



dt e 



iuit 



h{k,k')Gn^r.Ak',t) 



xGn^ny (k - k', t) + h2Gpp{k\ t)Gpp{k - k', t)\ (3.54) 

where Gpp{k,t), Gn^na:{k,t), Gnyny{k,t) are the correlation functions of the 
density and director modes respectively as functions of k and t. They 
can be found by inverse Laplace transforming the response functions. The 
momentum-dependent functions fo, fi, f 2, 13,90,9,92 /^land /12 are discussed 
in the Appendix ; their precise form is not necessary here. However, we 
do note that their angular dependence implies that all of the viscosities 
z/j,-i = 1,...5 are subject to density and director feedback. On the other 
hand, there is no feedback for c. A, and most importantly K/'ji. As we shall 
see in the next section, this latter result is important in eliminating the 
possibility of director freezing. 

4 Implications for Supercooling 

We now consider the implications of our results from the previous sec- 
tion for the supercooling of a nematic liquid crystal. Our analysis follows 
Leutheusser's original approach to the density feedback mechanism, except 
we are also interested in a potential director feedback mechanism, and the 
behavior of the six viscosity coefficients. 



1^ 



We begin by recalling Leutheusser's argument for the freezing of the den- 
sity fluctuations in a simple fluid. The density response function (3.11) can 
be rewritten as, 

^i{k,uj) = < 5p{k,uj)5p{~k,~Lj) > 

ujp+ik^r 

If the viscosity F grows as the fluid is supercooled (or compressed), then (4.1) 
indicates that 

^lik,^)^ 1 r^oo (4.2) 

UJ + tpC^/L 

This form indicates that the liquid is freezing, in particular as F ^ oo, $i de- 
velops a pole at cj = 0. In deriving (4.2) we assumed that F was growing very 
large. Leutheusser showed that this will in fact occur via a feedback mech- 
anism that couples density fluctuations to F, specifically the result (3.53), 
(dropping the terms proportional to G„„ for the moment). The correlation 
function Gpp{k,t) is given by the inverse Laplace transform of $1 multiplied 
by X- Thus, ignoring the k dependence, (3.53) yields on equation of the form, 

/•oo 

T{u) = F° + A / dt e*^*<I>?(t) (4.3) 

Jo 

where $i(t) is the inverse Laplace transform of $1(0;). Equations (4.2) and 
(4.4) can be solved to yield a glass transition where $1(0;) ~ l/cu, Gpp(t) -^ 
nonzero constant, as t ^ 00 and T{uj = 0) diverges. In particular F(c<j = 0) 
diverges as, 

T{lu = 0) ~ (T - Tor (4.4) 

where /i ~ 1.8 and Tq is the glass transition temperature. 

As discussed in the Introduction it is now believed that (4.2) is not correct 
and should be replaced by, 

^iik,uj)r^ I (4.5) 

where 7 does not go to zero as the fiuid is supercooled. The origin of 7 is 
still a subject of debate.®'^ Its presence will cutoff the Leutheusser transition 
and (4.4) will not be true asymptotically. Nevertheless, it is believed that 
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there will be a substantial growth in the viscosity as the fluid is supercooled 
which is eventually rounded off. 

We now proceed with a Leutheusser style analysis for the nematic, bearing 
the preceding discussion in mind about the limits of such an analysis. We 
first examine the longitudinal director response function given in (3.12) 



$2(^,^1;) = <5nx{k,uj)5nx{—k,—u> 

1 






(4.6) 



Lopo+ik'^yL 

where we have used (3.17). Equation (3.54) indicates that u^ does incorpo- 
rate density feedback (through the product Gpp{k' , t)Gpp{k — k', t)) and hence 
it will grow as the nematic is supercooled. Eqn. (4.6) reduces then in the 
large ul limit to: 

In the Appendix we show that K/'-fi does not renormalize in any dramatic 
way and is not affected by the density or possible director feedback mech- 
anism. Experimentally, -ft'/ 71 has been measured in supercooled nematics^° 
and found to be consistent with K/'ji ~ e~^/-^. Thus, as the temperature is 
reduced, the director mode will slow down but there will be no sharp transi- 
tion a ia Leutheusser or even a rounded one. The viscosity u^ will grow and 
appear to diverge due to the density fluctuations as will z/^. The absence of 
a freezing transition for the director modes is true even if the director modes 
are propagating rather than diffusive. In the limit of large ul, (4.7) is still 
obtained even if F^ and Tf in (3.17) become complex (which implies that 
there are damped, propagating modes). Thus to study the growth of the 
viscosities in (3.53) - (3.55), all terms proportional to either Gn^n^ or GnyUy 
can be dropped, and feedback from density fluctuations alone occurs. 

We now discuss the experimental implications of our results. First we 
summarize our results in terms of the viscosity coefficients z/i, z/2, z^s, z/4, z^s, 71 
and 72. The latter two viscosities as we have already indicated show no 
pretransitional growth, at least within the context of our theory. As noted 
above, the transport coefficient K/'-ji has been measured in supercooled ne- 
matics and only activated behavior is observed. However, we do expect that 
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the five viscosities z/j will show significant growth as the nematic is super- 
cooled. As they are all driven by the same density feedback mechanism we 
expect them to show the same power-law behavior indicated in (4.4), with 
eventual rounding off. At this time we cannot predict how wide the temper- 
ature regime will be where (4.4) is valid. We also expect on the basis of the 
results of refs. 18 and 21 where the effects of local structure on a simple fluid 
were considered that the shear modes of the nematic will become propagat- 
ing modes at sufficiently high frequency. If the Leutheusser scenario were 
not ultimately invalidated by a cutoff mechanism then at the glass transition 
and below the shear modes would propagate at any frequency and a nonzero 
shear modulus would be present. Again at this time we cannot predict what 
the lower frequency cutoff will be; calculation of this cutoff will be sensible 
once the controversy regarding the rounding off of the glass transition is re- 
solved. At that time it will also be sensible to study in detail (3.53) - (3.55) 
and calculate the k and u dependence of the viscosities. The modes associ- 
ated with the director relaxation remain diffusive. In that sense there will 
be an interchange of the "fast" and "slow" modes of eqns. (3.19) - (3.21), 
with director relaxation becoming fast and shear relaxation becoming slow. 

Table I. The zeroth order matrix Gq . The tensors aij and Lij are defined 
in eqns. (3.6) and (3.7), and x~^ = A + Bk^ 



Pi rij p Pj hj 



P 

Hi 
P 

hi 



-cj poX"^kj 

ki -uj6ij + iLij/po -j^Kaij 

-a,,e 5,,(^ + ^) 

uj -ki 

-PoX'^kj uSij + iLij/po ajik^ ^P^^Lij 

^a., S,,{u; + i^) 2/3" V' 
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Appendix 

In this appendix we provide the technical details of our one-loop calcula- 
tion of the self-energies appearing in (3.24) - (3.31). 

We consider only those diagrams where the propagators are either 
Gpp{k,uj) or Gnnik,uj) in order to search for a growth in the viscosities due 
to density and/or director feedback. We will ignore mixed propagators Gpn^ . 
As discussed by Forster, this correlation function contains both sound poles 
as well as Ts and T f. When expanded in terms of these poles, the strengths of 
the sound poles can be determined and for large viscosities yield a constant 
contribution, rather than the 1/uj needed for feedback. The strengths of the 
director poles cannot be determined to the order in k that we have worked, 
but it is unlikely that a 1/a; pole would emerge. 

The diagrams potentially yielding feedback are all of the bubble variety 
formed from two three-point vertices. In equations (3.42) - (3.52) we have 
listed those vertices which contribute to the self-energies we need. While we 
have chosen the external momentum k to lie in the x — y plane, we note that 
internal momenta can point in any direction. 

We find the following results for the self-energies with two hatted external 
lines: 



^P.P. = ^^{Gppik\.')Gp,ik-k',uj-.') 



(27r 27r 

■[^{Ak^ + B{k',{k-k') - kj' ■ {k-k'))) ■ 

{2Ak^B{k'^k'^ + kj' ■ ik' - k)) + 

^-f [kj. — k^)[kz — kj y[kx — k^)[kz — k^) + k^k^ 

-i\\ + i){^^{\ + i)klk'Ak.-k';) 

+kz\kx — k^)[kz — k^) ) • \k^[kz — k^) + k^[kx — k^)j 

-{liikx-k'jikz-k'^r+k'x') 

+ i{X + l)kz{kUkz-k',) + k',{kx-k'^))^ 

■{Akx + B{k'^{k ■ k') - kj' ■ {k - k')) 

+ Gn:,nAk', ^')Gn^n^ {k - k' , UJ - Lo') 
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■ {-K^Kik - k'fiKik - k'f + k'\k^ - t,)) 
~{2\K^k, + \'K^kl){k-k'f{k'' + {k-k'f} 

+Gnyny {k', Uj)Gnyny {k - k' , U - u') 

i-KXik-k'nk'Ak-kr + k'\k^ - k'^)} 



Pv"y 



dX duj' 



[G,p{k',uj')G,,{k-k\uj-uj') 



-^k'\k'^f{k.k') 



(27r)3 27r 
4 

+iX{k,-k':)\k'^{k':)'-k'^{k,-k':)') 

-/^(A + l)Q(A + l)A:X(A:.-A:;) 

— kzky{kz ~ kj j ■ {ky{kz — kj — kyk^j 

— \J[kyk^ — ky[kz — kj ) 

+L{x + i)k,{k'{k, - ki) - k'X)) 



■Bk'yik-k'] 

+Gn^n^ {k', Uj')Gn,n^ {k - k' , UJ - Uj) 

i-K^k'y{k-k'f{k'y{k-k'f-kX 

+ {Gnyny {k\ Uj)Gnyny {k - k' , UJ - Uj') 

■ (K^k'y{k -k'Y{k'y{k -k'f - k'yk 



''PzPz 



dX duj' 



|G'pp(A;', uj')Gpp{k -k',uj - uj') 



(27r)3 2tt 

-{Ak, + B{k',{k ■ k') - k,k' ■ {k - k)) 



■{2Ak, + B{k'X + kM' ■ (k' - k)) 
-l\k,-k'Mk.-k':)Xk'!) 
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-liih - Kf) + k'!) 

■{Ak, + B{k'^{k - k') - kS' ■ {k - kf)) 

+ [Gn^nAk\^)Gn^nA^ - k' ,UJ - Uj') + 
Gnyny{k'uj')Gnyny{k - k' ,UJ - Cj')) 

i-K'){kAk-kr{K{k-k)' 
+ k'\k-,-k'S)+{k'^ + {k-k'f). 
{k-k'f{\^kl-2\k,k'S) 



'p^p. 



I JW^ {Gppik', Lo')Gpp{k -k',uj- Lo') 
[Ak^ + B{k'Ak ■ k') - kj' ■ {k - k'))) 
■ (2Ak, + B{k',k'^ + k,k' ■ (k' - k))) 

-LiAk. + Bmk.k')-kJ'.Ck-k')) 
{{k, - k'^f + kf) + (A - l)kX{K - ki)) 

-i\k. - k'^){K - k':)\{K - k':)' + kT 



-(A + i)kAk. - k';)\k'Ak. - k'^) + k'Ak. - k',)) 



Lf 



l\\-l)k,kAh,-k',){K-k 

p 

— — (A — l)kzkzk^k^[kz — k^) 

-L[\ + i)kzkAk^-k',). 

{Akz + B{k'Ak ■ k') - kzk' ■ {k - k') 

+Gn^n^ {k', Uj')Gn,n^ - k' , UJ - Uj') 

■{k - k'fK\-]^{\ + l)kAk'z{k - k')' + k'\kz - k'^)) 
-\{\-l)kMk" + {k-k'f)) 
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+\^K^k^k,{k-k'f{k'^ + {k-k'f) 

+ {Gn^nAk',i^')Gn^nAk - k' ,UJ - Uj')) 
+Gnyny {k', Uj')Gnyny {k - k' , UJ - Uj')) 

.{-K'){k-k')X 
(k'^{k-k'f + k'\k,-k',) 

+k,k'^ + ik- k'fk,) 



Sn^ft. = j j^-^[Gppik\uj')Gpp{k-k\uj-J) 

■( ) k^[kz — kj{k^{k^ — kz) + k^[k^ — kx)) 

71 

The remaining self-energies needed, S^^j^j;, '^■hyPy and Sp , do not exhibit 
feedback. The first of these renormahze the Frank constants and has no 
graphs of the form we are considering. However, the finite graph contributing 
to T^n^fix does break the one-constant approximation. Similarly, the other two 
self-energies also do not have any graphical contributions of the form we are 
considering. 

With the self-energies (Al) - (A5) the viscosities Viii = 1, ... 5 can prin- 
ciple be calculated to one-loop order using (3.32) - (3.36). In general these 
expressions are quite complicated and not much will be learned by display- 
ing them. Even the hydrodynamic limit is difficult to evaluate (except for 
V2 and 2/3) due to the anistropy of the propagators. Finally the functions 
/o, /i, /2, /s, fi'o, 5'i, fi'2, hi and /12 appearing in (3.53) - (3.55) can in principle 
be calculated using (3.37) - (3.39) and (Al) - (A5). 
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